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ABSTRACT 

We combine surveys of the radio sky at frequencies 22 MHz to 1.4 GHz with data from the ARCADE- 
2 instrument at frequencies 3 to 10 GHz to characterize the frequency spectrum of diffuse synchrotron 
emission in the Galaxy. The radio spectrum steepens with frequency from 22 MHz to 10 GHz. The 
projected spectral index at 23 GHz derived from the low-frequency data agrees well with independent 
measurements using only data at frequencies 23 GHz and above. Comparing the spectral index at 
23 GHz to the value from previously published analyses allows extension of the model to higher 
frequencies. The combined data are consistent with a power-law index f3 = —2. 64 ±0.03 at 0.31 
GHz, steepening by an amount A/3 = 0.07 every octave in frequency. Comparison of the radio data to 
models including the cosmic ray energy spectrum suggests that any break in the synchrotron spectrum 
must occur at frequencies above 23 GHz. 

Subject headings: radio continuum: general, radiation mechanisms: non-thermal 



1. INTRODUCTION 

Synchrotron emission from relativistic cosmic ray elec- 
trons accelerated in the Galactic magnetic field domi- 
nates the diffuse radio continuum at frequencies below 
1 GHz. It is an important foreground contaminant for 
measurements of the cosmic microwave background ra- 
diation, and also serves to probe the Galactic magnetic 
field and cosmic ray distributions. Measurements of the 
synchrotron frequency spectrum are thus of interest to 
several areas in astrophysics. 

An isotropic distribution of relativistic electrons at a 
single energy E = jmc 2 propagating in a uniform mag- 
netic field B has emissivity 
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where a is the pitch angle between the magnetic field 
and the line of sight, and 



F{x) = x K b/ z(x')dx' 
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is defined in terms of the modified Bessel function of 
order 5/3 with variable x = vjv c and 
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(jSchwinger I fl94l iWestfold et aT~l fl95l IQster I fl96ll) . 
For a power-law distribution of electron energies 
N(E) oc E p propagating in a uniform magnetic field, 
the synchrotron emission is also a power law, 
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where Ta is antenna temperature, v is the radiation fre- 
quency, and 

P= P -^ (5) 
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(jRvbicki fc Lightman Ill979f ). 

Measurements of the synchrotron spectral index pro- 
vide important input for models of cosmic ray propa- 
gation. Solar modulation reduces the local cosmic ray 
electron density for electron energies below a few GeV so 
that synchrotron emission provides the most direct probe 
of low-energy cosmic rays. Measurements of the cosmic 
ray spectrum above a few GeV in turn inform models of 
the high-frequency synchrotron spectrum. Energy losses 
from cosmic ray propagation steepen the cosmic ray spec- 
trum, increasing p toward higher energies. The observed 
steepening from p ~ —2.6 at 5 GeV to p ~ —3.2 at 
50 GeV predicts a corresponding steepening in the syn- 
chrotron spectru m from (3 ~ —2.8 at 1 GHz to P ~ —3.1 
at 100 GHz (jStrong. Moskalenko. fc Ptuskin ll2007h . 

Comparison of the cosmic ray spectra to the predicted 
synchrotron spectrum is complicated by confusion from 
competing radio emission sources. The diffuse radio con- 
tinuum is a superposition of the cosmic microwave back- 
ground, synchrotron emission, free- free emission from the 
warm ionized interstellar medium, and emission from in- 
terstellar dust. A number of authors have attempted to 
disentangle the various emission sources to determine the 
synchrotron spectral index (for a recent revi ew see Ap- 
pendix A of iStrong, Orlando, fc Jaffe I (|2011[ )). Despite 
some discrepant results, the general trend shows a steep- 
ening of the synchrotron spectrum from f3 ~ —2.5 at 22 
MHz to f3 ~ —3.0 above 23 GHz, in rough agreement 
with the observed cosmic ray spectra. 

Several factors contribute to the observed scatter in 
estimates of the synchrotron spectral index. Most es- 
timates, particularly those below 23 GHz, assume a 
power-law spectrum for synchrotron and do not explic- 
itly model spectral steepening. Comparisons between 
closely-separated frequencies more accurately reflect the 
local synchrotron spectrum, but have larger uncertain- 
ties from competing emission sources or measurement 
offsets. Analyses with broader frequency coverage reduce 
foreground and offset uncertainties but average over any 
spectral steepening. 

Two additional effects are important for analy- 
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ses including data from the Wilkinson Microwave 
Anisotropy Probe (WMAP) at frequencies 23 to 94 
GHz. A growing body of evidence suggests that a 
substantial fraction of the diffuse continuum near 
23 GHz consists of electric dipole radiation from a 
population of small, rapidly spinning; dust grains 
JlKogut et al 1996; de Oliveira-Costa et al. II1997L 12004 
Miville-Deschenes et al. I 120081: Dobler & Finkb einer I 



2008; 



2010; 



lYsard. Miville-Deschenes. fc Verstraete I 

iKogut et al I [20111 : IGold et al. I EnTF 
Planck co llaborat ion I l2011f ). Spinning dust emis- 
sion is expected to peak at frequencies near 2 GHz 
dDraine fc Lazarian I fl998L PAli-Hai'm oiid et al. I 12001 
Hoang. Draine. fc Lazarian ] [201(1 [Ysard & Verstraete 



201 Of ) . Analyses that ignore this component to attribute 



the observed emission only to synchrotron radiation tend 
to over-predict the synchrotron amplitude at frequencies 
near 23 GHz, biasing the derived spectral index to 
flatter values when comparing to lower frequencies and 
steeper values when comparing to higher frequencies. 

A second systematic error can result from improper 
treatment of offsets in the data. Measurements from ra- 
dio surveys at frequencies below 20 GHz generally include 
the absolute intensity (zero level) of the sky. The WMAP 
differential radiometers are insensitive to any constant 
(monopole) intensity on the sky; the zero level of the 
WMAP sky maps is set so that the map intensity in the 
Galactic pol ar caps matches a cosecant fit to the mid - 
latitude sky (jBennett et al. 1120031 : iHinshaw et al. II2009D . 
Analyses that directly compare low-frequency radio sur- 
veys to the WMAP data without subtracting a monopole 
component from the radio data will miss the correspond- 
ing emission in the WMAP bands, biasing the derived 
spectral index to steeper values. 

The zero l e vel b ias can be significant. The 
Hasl am et al. I (pj)81) survey at 408 MHz is commonly 
used to model synchrotron emission. The North Galac- 
tic pole has measured temperature 19 ±3 K at 408 MHz, 
while a esc \ b\ fit to the same 408 MHz map predicts a 
polar contribution of only 5.1 ± 0.6 K. Similar results 
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Fig. 1. — Toy model showing effect of improper zero level sub- 
traction on the derived spectral index between 408 MHz and 23 
GHz. 



TABLE 1 

Sky Surveys Used for Synchrotron Analysis 

Frequency Calibration Offset Relative 

(GHz) Uncertainty Uncertainty (K) Uncertainty 61 
O022 0^5 5000 0T5 

0.045 0.10 250 0.11 

0.408 0.10 3.0 0.17 

1.420 0.05 0.5 0.63 

3.20 0.001 0.011 0.10 

3.41 0.001 0.006 0.07 

7.98 0.001 0.036 0.89 

8.33 0.001 0.042 2.64 

9.72 0.001 0.003 0.34 

10.49 0.001 0.002 0.27 

a Quadrature sum of calibration and offset uncertainties, 
divided by the mid-latitude sky temperature. 

apply to the Southern hemisphere, where the measured 
polar cap temperature of 21 ± 3 K significantly exceeds 
the value 4.0 ±0.5 K obtained from a esc |6| fit. Only 2.7 
K of the difference can be attributed to emission from the 
cosmic microwave background, leaving a large residual. 

Figure [T] illustrates the bias induced by including this 
residual at 408 MHz but excluding it from the WMAP 
data. We take the 408 MHz map, remove the 2.7 K 
CMB monopole, and scale the remaining (radio) emis- 
sion to 23 GHz using a power-law index j3 = —2.7 to 
produce a toy model of the radio sky at 23 GHz. Follow- 
ing the WMAP processing, we then remove a monopole 
from the scaled map so that the map temperature in the 
south polar cap matches the esc |6| fit. We then com- 
pute the bias in apparent spectral index by comparing 
the 408 MHz map to the scaled 23 GHz map before and 
after removing the scaled monopole. Figure [I] shows the 
bias in spectral index at 23 GHz, binned by Galactic lati- 
tude. Dis-similar treatment of the map zero level creates 
a spatially varying bias A/3 ~ 0.15, comparable to the 
total spectral steepening predicted by the measured cos- 
mic ray spectra. The bias is largest in regions where the 
sky brightness is faintest, at high latitudes or away from 
the Galactic center. 

Measurement uncertainties in the absolute level of the 
sky brightness can also introduce bias in estimates of the 
synchrotron spectral index. Many of the low-frequency 
surveys have uncertainty in the measured zero level ap- 
proaching 30% of the polar cap brightness. As with 
the toy model above, such measurement errors intro- 
duce spatially dependent biases that are largest where 
the sky brightness is faintest. Minimizing uncertainty in 
the derived synchrotron spectrum requires a combination 
of measurements with good sky coverage and good con- 
trol of offset uncertainty at frequencies where competing 
emission sources are faint. No such ideal data set yet 
exists. In this paper, we model the synchrotron spectral 
index and curvature using low-frequency radio surveys 
with high sky coverage but large offset uncertainty, com- 
bined with higher-frequency measurements with limited 
sky coverage but still useful offset uncertainty. 

2. SKY MAPS 

We mode l synchrotron em ission using radio surveys 
at 22 MHz (iRoger et al. 111999ft . 45 MHz dMaeda et al ' 
19991: lAlvarez et al. I 119971), 408 MHz ([Haslam et al 
1981 ), and 142 MHz (jReich. Testori. fc Reich I 120011 : 
Reich & Reich 1986). These surveys have full or nearly- 
full sky coverage at frequencies where Galactic radio 
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Fig. 2. — Sky coverage for this analysis. The plot shows the 408 
MHz sky survey convolved to 11? 6 angular resolution. Pixels com- 
mon to all 10 radio surveys are shown in color. The sky coverage is 
limited by the ARCADE 2 observations but includes the Galactic 
plane, mid-latitude sky, and portions of the North Galactic Spur 
(radio Loop I). 

emission is significant, with gain and zero- level systemat- 
ics controlled at the 10-20% level. We supplement the ra- 
dio surveys with sky maps from the Absolute Radiometer 
for Cosmology, Astrophysics, and Diffuse Emission (AR- 
CADE 2) instrument^ at 3, 8, and 10 GHz (jKogut et al. I 
120111 1. The ARCADE 2 data observe both the Galac- 
tic plane and mid-latitude regions (|6| < 40°) with suf- 
ficient control of zero-level uncertainty to constrain the 
synchrotron curvature relative to the lower-frequency ra- 
dio surveys. 

Table 1 summarizes the input sky maps. The increase 
in the offset uncertainty at low frequency is compensated 
by a corresponding increase in sky brightness. The fi- 
nal column shows the relative measurement uncertainty 
for a mid-latitude region, defined as the ratio of the 
combined offset and calibration uncertainty to the mea- 
sured brightness at (/, b) = (17°,— 35°) after removing 
the CMB monopole. The selected maps provide roughly 
uniform relative sensitivity to synchrotron emission over 
2.5 decades of frequency. 

We convert all maps to units of antenna temperature 
and subtract the CMB monopole at (thermodynamic) 
temperature 2.725 K from the measured sky tempera- 
tures. We then convolve each map to the 11? 6 angular 
resolution of the ARCADE 2 instrument. At frequen- 
cies of 10 GHz and below, both thermal dust emission 
and spinning dust emission are negligible. Free- free emis- 
sion, however, can still be appreciable. We correct the 
convolved maps by scalin g the WMAP 7- ye ar maximum 
entropy free-free model (Gok Tetal. Il201l| ) to each fre- 
quency using spectral index —2.15, convolving the scaled 
model to 11?6 angular resolution, and subtracting the 
convolved model from each sky survey. The resulting 
maps are dominated by synchrotron emission. 

3. ANALYSIS 

2 The ARCADE data are available at the Legacy 
Archi ve for Microwave Background Data Analysis, 
http : / /lambda . gsf c . nasa . gov 



The input sky maps define a data set T(n, v) sampled 
at discrete pixel directions n and 10 discrete frequencies 
v ranging from 22 MHz to 10 GHz. We restrict the anal- 
ysis to the 8% of the sky observed at all 10 frequencies. 
Figure [2] shows the resulting sky coverage. Within the 
common sky coverage, we model synchrotron emission as 
a modified power law 



T{n. v) = A(n) ( ^ 



(6) 



with spectral index j3 and curvature C defined with re- 
spect to reference frequency v$ = 310 MHz. The adopted 
value for minimizes covariance between the fitted am- 
plitude A and spectral index /3, simplifying extrapolation 
to other frequencies. 

For each pixel h we define a 10 by 10 data covariance 
matrix M with diagonal elements determined by the in- 
strument noise, calibration and offset uncertainty. The 
input maps are not all linearly independent. Measure- 
ments at 22 MHz used the 408 MHz map to determine 
the declination dependence of the gain. We estimate the 
resulting correlation of spatial structure in the two maps 
at 50%. The ARCADE 2 maps have independent instru- 
ment noise but share a fraction of the offset uncertainty 
relate d to absolute therm ometry uncertainty and ground 
glint (jSingal et al. Il20llh . All 10 maps share a common 
model for free-free emission. We conservatively estimate 
the uncertainty in the free-free correction at 30% of the 
free- free amplitude; however, the results do not change 
significantly as the model free-free amplitude is varied 
by as much as 50%. Off-diagonal elements in M include 
these effects. 

For each pixel, a least-squares minimization determines 
the best-fit parameters A, /3, and C. Figure [3] shows the 
measured temperature and best-fit model for the bright- 
est Galactic plane region (Z,6) = (52°, 0°) within the 
common sky coverage. The data show evidence for spec- 
tral curvature, with best-fit values j3 = —2.60 ±0.04 and 
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Fig. 3. — Sky temperatures and best-fit model (solid line) for a 
4° diameter patch on the Galactic plane centered on the the bright- 
est pixels in the ARCADE 2 sky coverage. For clarity, the data are 
plotted relative to the best-fit power-law model ((3 = —2.56) with 
zero curvature (dashed line). All fits include the non-trivial covari- 
ance between individual data points. The data are well described 
by a model with spectral index (5 = —2.60 ± 0.04 and curvature 
C = -0.081 ±0.028. 
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Fig. 4. — Sky maps of best-fit spectral parameters. Left to right: synchrotron amplitude, spectral index, and curvature evaluated at 
reference frequency 310 MHz. The top panels show results from the 10-frequency fit, while the bottom panels include the spectral constraint 
at 23 GHz. 

C = -0.081 ± 0.028 evaluated at v = 310 MHz. The 
spectral curvature C for this region is significant at ap- 
proximately 3 standard deviations compared to the base- 
line model with (7 = 0. 

We may extrapolate the spectral models to compare 
the results at frequencies 10 GHz and below to indepen- 
dent determinations of the spectral index using WMAP 
data at frequencies 23 GHz and above. We compute 
the antenna temperature of the modeled spectra to de- 
rive the effective power-law index for frequencies near 23 
GHz. Note that this is not equivalent to evaluating Eq. 
[6] at v — 23 GHz, which would yield the scaling from 310 
MHz to 23 GHz but not the power-law index at 23 GHz. 
The mean for all 258 pixels in the common sky coverage 
is /#23 = —3.02 with standard deviation 0.22. 

The extrapolated value compares well with indepen- 
dent determinations of the spectral index above 23 GHz. 
iKogut et all (|2004f ) analyze WMAP polarization data 
to derive synchrotron spect ral index /3 = —3.2 db 0.1 
averaged over the full sky. iDunklev et all (2009) use 
a Bayesian analysis of polarization data and find the 
mean spectral index j3 = —3.03 ± 0.04 with pixel-to- 



pixel standard dev iation 0.25 over the high-latitude sky. 
iGold et al. I (|2011[ ) use template fitting techniques to de- 
rive spectral index f3 = —3.13 between 23 and 33 GHz. 
The mean spectral index derived from the 10 low- 



frequency radio surveys agrees with the value derived 
from WMAP data at higher frequencies. Much of the 
scatter in the extrapolated spectral indices results from 
pixels at high latitude where the emission is faintest. 
The extrapolated index in these pixels can reach un- 
physical values. We reduce the scatter in the fitted 
spectral parameters by applying additional constraints 
using WMAP data at 23 GHz and above. The sim- 
plest such constraint, adding WMAP temperature data 
to the multi- frequency fit, is problematic. Not only would 
the procedure need to include additional free parame- 
ters to account for emission from thermal dust or spin- 
ning dust (both negligible at lower frequencies), but each 
of the low-frequency maps would require a correction 
to remove the monopole contribution missing from the 
WMAP data. Although the WMAP zero level is clearly 
defined by a esc \ b\ fit to mid-latitude data, the astro- 
physical interpretation of a similar procedure applied to 
low- frequency radio surveys is less clear. The coldest re- 
gion of the radio sky is not at the Galactic poles, but at 
mid-latitudes above the Galactic anti-center. Subtrac- 
tion of too large a monopole can leave unphysical neg- 
ative residuals. Limited sky coverage exacerbates this 
problem. 

We avoid these problems by using a constraint based 
on the spectral index derived solely from WMAP data. 
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For each pixel, we use the radio data (Table [T]) to fit the 
synchrotron amplitude A(h) over a 2-dimensional grid in 
the spectral parameters $ and C (Eq. [6j). At each grid 
point, we compute the x 2 value R T M~ 1 R where M _1 
is the inverse covariance matrix and R is the difference 
vector between the measured and modeled temperatures. 
We then use the spectral parameters /3 and C to evaluate 
the power-law index at 23 GHz and compare the result- 
ing value to a prior. We use the difference between the 
extrapolated spectral index and the prior to augment the 
X 2 at each grid point, 



2 2 

x 



(7) 



where ^23 is the model spectral index evaluated at 23 
GHz, and j3 p ± a v — — 3.1 ± 0.1 is the prior at 23 GHz. 
The minimum x 2 over the entire grid then defines the 
best-fit model at that pixel. This allows inclusion of the 
spectral information derived from frequencies above 23 
GHz without confusion from either additional emission 
components (thermal or spinning dust) above 23 GHz or 
the missing zero level in the WMAP data. 

Figure |4] shows the best-fit spectral parameter maps, 
while Figures [5] and [6] show the distribution of the best- 
fit values for the spectral index and curvature. Including 
the constraint at 23 GHz, the best-fit spectral index has 
mean value —2.64 and standard deviation 0.03 at refer- 
ence frequency = 310 MHz. The best-fit curvature has 
mean value —0.052 with standard deviation 0.005. The 
corresponding spectral index at 23 GHz is (^23) = —3.09 
with standard deviation 0.05. Inclusion of the prior at 
23 GHz does not induce a significant shift in the mean 
for the extrapolated spectral index, but does significantly 
reduce the pixel-to-pixel scatter. 

Comparison of the parameter distributions with and 
without the spectral constraint at 23 GHz demonstrates 
that the addition of the spectral constraint mainly affects 
the fitted curvature values (Figs. [5] and [6]). The adopted 
value for the spectral constraint comes from independent 
analysis of WMAP data and is weighted toward regions 
of higher synchrotron intensity. We test whether this 
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Fig. 5. — Spectral index evaluated at 310 MHz. The dashed line 
shows the distribution from the 10-frequency radio data while the 
solid line includes the prior at 23 GHz. 
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Fig. 6. — Spectral curvature at 310 MHz. The dashed line shows 
the distribution from the 10-frequency radio data while the solid 
line includes the prior at 23 GHz. The mean curvature C = —0.052 
corresponds to a steepening of the local spectral index by an 
amount A/3 = 0.07 every octave in frequency. 



creates a bias in the fitted curvature values by splitting 
the observed sky coverage into two subsets of equal area, 
defined by the brightest and faintest 50% of the fitted 
amplitudes at 310 MHz. Within each subset, we com- 
pare the mean and standard deviation of the fitted cur- 
vature derived from the 10-frequency fit or the enhanced 
fit including the constraint at 23 GHz. The "bright" sub- 
set shows no shift in the mean spectral curvature when 
the 23 GHz constraint is added (although the scatter 
is significantly reduced). The 'faint" subset shows both 
a reduction in scatter and a modest shift in the mean 
value, with curvature parameter steepening from -0.034 
to -0.049 when the 23 GHz constraint is included. This 
shift is less than one standard deviation: given the lim- 
ited sky coverage, the available radio data do not yet 
provide significant evidence for spatial variation in the 
synchrotron curvature. 

4. DISCUSSION 

Radio data show statistically significant steepening of 
the synchrotron spectrum from 22 MHz to 10 GHz. The 
nearly uniform relative uncertainty of the selected data 
minimizes dependence of the fitted parameters on off- 
set or calibration errors at any one frequency. We test 
whether the best-fit parameters are particularly sensi- 
tive to any one input map by repeating the analysis af- 
ter dropping one or two maps from the fit. We select 
either one radio survey (22 MHz, 45 MHz, 408 MHz, 
or 1420 MHz) or a pair of ARCADE frequency chan- 
nels (3 GHz, 8 GHz, or 10 GHz) and repeat the fit af- 
ter deleting the corresponding elements from the data 
vector T and covariance matrix M. The resulting shift 
in either the spectral index or curvature parameters is 
smaller than the pixel-to-pixel standard deviation using 
all 10 frequency channels. Systematic errors in the offset 
or temperature calibration do not appear to dominate 
the multi- frequency analysis. 

The steepening of the synchrotron spectrum is broadly 
consistent with models of cosmic ray propa gation in the 
Galactic magnetic field. iJaffe et al. I (|2Qllh combine ra- 
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TABLE 2 
Local Power-Law Spectral Index 



Frequency 


Power- Law 


(GHz) 


Index 


0.022 


-2.36 


0.045 


-2.44 


0.408 


-2.67 


3.3 


-2.89 


23 


-3.09 


33 


-3.13 


41 


-3.15 


61 


-3.19 


94 


-3.24 



dio observations with the GALPROF0 cosmic ray propaga- 
tion code to model synchrotron emission on the Galactic 
plane. They find a power-law index —2.8 < f3 < —2.74 
from 408 MHz to 2.3 GHz and -2.98 < /3 < -2.91 from 
2.3 GHz to 23 GHz. The corresponding values for the 
spectral steepening model presented here are j3 = —2.76 
from 408 MHz to 2.3 GHz and fi = -2.97 from 2.3 GHz 
to 23 GHz, in good agreement with the cosmic ray model. 

We may use the best-fit values in each pixel to predict 
the synchrotron spectrum at higher frequencies where the 
emission is fainter and competing sources stronger. Pre- 
vious attempts to disentangle competing emission from 
free- free, synchrotron, thermal dust, and spinning dust 
emission using WMAP data have suffered from degener- 
acy between the synchrotron and spinning dust emission, 
both of which are falling at frequencies abo ve 33 GHz 
(see, e.g., the discussion in lGold et al. I (j2011f )). Extend- 
ing the synchrotron curvature observed at lower frequen- 
cies into the millimeter band reduces confusion between 
the spinning dust and synchrotron spectra and may fa- 
cilitate characterization of both the spatial distribution 
and frequency spectrum of spinning dust emission in the 
interstellar medium. 

Table [2] shows the local power-law index T oc v& at se- 
lected frequency bands. The modeled spectrum steepens 
by A/3 = 0.07 every octave in frequency, from f3 = —2.67 
at 408 MHz to /3 = -3.24 at 94 GHz. Note, however, 
that the spectral steepening observed at low frequencies 
can not continue indefinitely. Fermi measurements of the 
cosmic ray energy spectrum are consistent with a single 
powe r law from energy 7 GeV to 1 TeV (|Abdo et al. I 
120091 : lAckermann et al. 1 120101 . I2Q12D . If anything, the 
Fermi data suggest a modest flattening of the cosmic 
ray energy spectrum at higher energies, which would in- 
duce a positive curvature to the synchrotron spectrum 
a t frequencies above 23 GHz . 

Str ong. Orlando, fc Jaffe I (|2011l ) combine radio data 
at 22, 45, 150, 408, and 1420 MHz with WMAP data 
at 23 through 94 GHz and Fermi Large Area Telescope 
cosmic ray measurements and the galprop code to esti- 
mate the magnetic field intensity and synchrotron spec- 
trum from 22 MHz to 94 GHz. Figure compares the re- 
sulting synchrotron spectral index (diffusion model with 
injection index 1.3) to the curvature model from this pa- 
per. The cosmic ray model analyzes a limited latitude 
range 10° < \b\ < 50°, which we follow in Figure by 
excluding pixels at latitudes \b\ < 10°. The two methods 
agree for frequencies below 408 MHz (where they share 
common radio data) but differ at higher frequencies. The 

3 http://galprop.stanford.edu 
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Fig. 7. — Comparison of the synchrotron spectral index as a 
function of frequency for different models. Th e break in the 
spectral index for the Strong, Orlando, & Jaffe (|201lD cosmic 
ray model at frequencies of a few GHz is not reproduced by the 
ARCADE 2 observations at 3-10 GHz, which prefer models with 
nearly constant spectral curvature. 



cosmic ray model shows an increased synchrotron curva- 
ture from 408 MHz to a few GHz, followed by a spectral 
break to near-constant index j3 ~ —3.1 at higher frequen- 
cies. The ARCADE 2 data at 3-10 GHz do not reproduce 
these features, but are instead consistent with constant 
spectral curvature from 22 MHz to 10 GHz (Fig ED- 

Both models agree at frequencies near 23 GHz. The 
spectral index at 23 GHz derived from the 10-frequency 
radio data without the 23 GHz prior is consistent with in- 
dependent measurements of the index above 23 GHz and 
with the full radio fit including the 23 GHz prior. The 
radio data, taken alone, do not support a break in the 
synchrotron spectrum at GHz frequencies. Comparison 
of the radio fit to the cosmic ray model suggests that any 
spectral break in the synchrotron spectrum must occur 
at frequencies above 23 GHz. Direct confirmation of the 
synchrotron spectrum above 23 GHz remains a challenge. 



5. CONCLUSIONS 

Radio data are consistent with a synchrotron spectrum 
that steepens with frequency from 22 MHz to 10 GHz. 
Direct comparison of low-frequency radio surveys with 
the WMAP data at 23 to 94 GHz is complicated both 
by the presence of additional emission components at 
higher frequencies and by the subtraction of a substan- 
tial monopole component of sky emission by the differen- 
tial WMAP instrument. The synchrotron spectral index 
at 23 GHz, derived using only lower- frequency radio sur- 
veys, is consistent with the value derived independently 
using only data at higher frequencies. We extend the ra- 
dio data by comparing the extrapolated index at 23 GHz 
to a prior based on higher-frequency data. The combined 
data have mean spectral index f3 = —2.64 ±0.03 and cur- 
vature C = —0.052 ± 0.005 at reference frequency 0.31 
GHz. The measured spectrum steepens by an amount 
A/3 = 0.07 every octave in frequency. Comparison of 
the radio data to models including the cosmic ray en- 
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ergy spectrum suggests that any break in the synchrotron 
spectrum must occur at frequencies above 23 GHz. 
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